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Measurements are recorded for multiwire thermocouples consisting of either two or three wires of 
unequal diameters. Signals from the multiwire probe are recorded for a reversing gas flow with both 
a periodic temperature and time constant fluctuation. It is demonstrated that the reconstructed signal 
from the multiwire thermocouple requires no compensation provided o>/&>,<2.3 for two wires or 
3.6 for three wires where o>] ( = 27 r/) is the natural frequency of the smaller wire based on 
the maximum gas velocity. The latter results were possible provided Fourier transformed data from 
the wires were used and knowledge of the gas velocity phase angle was available. © 1995 
American Institute of Physics. 


I. INTRODUCTION 

The advantages of thermocouples are their low cost, re- 
liability, and simplicity since they do not require optical ac- 
cess or elaborate support electronics. The design of a ther- 
mocouple represents a compromise between accuracy, 
ruggedness, and rapidity of response. Although thermo- 
couples are suitable for the measurement of high-frequency 
temperature fluctuations (<1 KHz) in a flowing gas or liq- 
uid, the measured signal must be compensated since the fre- 
quency of the time-dependent fluid temperature can be much 
higher than the natural frequency of the thermocouple 
probe. 1 Moreover, use of a single wire thermocouple requires 
knowledge of the fluid velocity (normally assumed constant) 
and fluid properties (e.g., viscosity, density, etc.) to deter- 
mine the convective heat transfer coefficient of the wire and 
its natural frequency. 

The present article describes the performance of novel 
two or three-wire thermocouples with unequal wire diam- 
eters for use in unsteady or irregular flows of varying veloc- 
ity. In this case, the time constant or natural frequency of the 
thermocouple wire is time dependent and the normal proce- 
dure for compensation of high-frequency signals is imprac- 
tical. Previous work 2-4 indicated the usefulness of the con- 
cept in constant velocity flows, where it was demonstrated 
that no compensation was required nor any knowledge of the 
fluid velocity or properties over a useful range of fluid tem- 
perature frequencies. In the present article, the multiwire 
thermocouple has been tested in a reversing flow field. The 
results of experimental measurements are presented along 
with the suggested procedure for the reduction of the data 
from the multiwire thermocouple, as shown in Fig. 1 . 


II. THEORY 


lar fluid motion are (a) cellular vorticies in stratified fluids or 
(b) oscillating flows in the exhaust from certain engines (e.g., 
Stirling, etc.). 

A. Natural frequency of single wire 

The conservation of energy for a single wire is 
dT 

~dt ==0) rx{T g -T), ( 1 ) 

where the natural frequency in radians/s 
Ah 

<S>n ~~pcD' 

Here: 

(o n = natural frequency (s~ 1 ) , 

T= temperature of the thermocouple(°K), 

T g = temperature of the gas(°K), 
t=time(s), 

p = wire density (kg/m 3 ), 
c = heat capacity of wire(J/kg— °K), 
h = heat transfer coefficient W/m 2 — °K), 

D — wire diameter(m). 

It is convenient to rewrite the heat transfer coefficient in 
terms of the Nusselt number 

kg 

Nu, (3) 

where 

kg=thermal conductivity of the gas(W/m— °K) 


Unsteady fluid motion develops in a variety of circum- 
stances. For example, surfaces of discontinuity behind sharp 
edges such as airfoils and adverse pressure gradients on blunt 
objects and in sharp corners lead to flow separation and the 
formation of recirculating eddies. 5 Other examples of irregu- 


Nu=Nusselt number. 

Thus, the natural frequency in Eq. (2) becomes 
_ 4/igNu 


(4) 
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FIG. 1. Schematic of a three- wire thermocouple. Wire diameters are 2, 4, 
and 8 mils. 

where 6 

Nu=C 0 Re^Pr 173 ( 5 ) 


TABLE I. Three-wire diameters in Fig. I. 

) D (m) £>(mil) 


52.6 

50.8X10" 6 

18.3 

101.6XI0“ 6 

6.4 

203.2X10 -6 


Substituting values corresponding to a type K thermocouple 
and a maximum gas velocity V max =25 m/s -1 for air at stan- 
dard conditions into Eqs. (7) and (8), one obtains values of 
the maximum natural frequency o> max . 


B. Two wires 

The conservation of energy given by Eq. (1) for a single 
wire contains two unknowns: o) n the natural frequency and 
the desired gas temperature T g . Each additional wire of a 
different diameter adds an unknown natural frequency co n . 
We also assume for a multiwire thermocouple that additional 
information is available concerning the ratios of natural fre- 
quencies. The latter provides a sufficient number of indepen- 
dent equations to determine the time-dependent gas tempera- 
ture T g . 

Consider two thermocouple wires of unequal diameter 
D\<D 2 . The conservation of energy for each wire becomes 

dT i 

- 5r = c|v|"Dr 2 (r g -r 1 ), (9) 


and . * 


Re= 


\V{t)\D 


Here, Re and Pr are the Reynolds and Prandtl numbers, re- 
spectively, and 

fi = viscosity of gas(kg/ms), 

v g = kinematic viscosity of gas(m 2 — s -1 ), 

V(?) — time dependent fluid velocity (m/s). 

The natural frequency of a single wire from Eqs. (4) and 
(5) can now be written in the form, 

k g C 0 4 Re m Pr ,/3 

“* w (6) 

or separating the wire diameter and fluid velocity 

u n = C\V\ m D m ~ 2 , (7) 

where the coefficient C depends on both fluid and wire pa- 
rameters or 


4*„C 0 Pr 1/3 

~{pc)v m g 


( 8 ) 


dT 2 

dt 


C\V\”DT 2 (T g ~T 1 ). 


( 10 ) 


In principle, knowledge of the wire diameters provides an 
additional equation 


( o ] _(D l \ m ~ 2 
W 2 \T>2 / 


( 11 ) 


representing the ratio of natural frequencies as defined by Eq. 
(7). Dividing Eq. (9) by Eq. (11) and solving for the gas 
temperature one obtains 



where T x =dT x /dt. 

Our experience with multiwire thermocouples in con- 
stant flow 3,4 and the present study indicate that construction 
of the gas temperature from Eq. (12), which requires the 
measured wire signals T x , T 2 and their derivatives, is un- 
stable and unreliable. We have found in all cases, however, 
that satisfactory results are obtained by taking the Fourier 
transform of the conservation equations and deriving the gas 
temperature in terms c0 transformed variables. For the case 
of unsteady flow, independent knowledge of the gas velocity 
is therefore required. 
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Thus, taking the fast Fourier transform (FFT) of Eqs. (9) 
and (10) and dividing the expressions, the transform of the 
product of gas temperature and fluid velocity is obtained in 
the form 


T = 

1 gv 




m — 2 


m - 2 


(13) 


T 2 j \D 2i 

where f 1 „=FFT(|V r |“7’ l ) and FFT (\V\ m T g ). For con- 
venience we also define the parameter /3=w 1 /u) 2 or 

i m — 2 

(14) 


I D , 


The inverse transform of T gv provides the reconstructed gas 
temperature or 


1 

Tg ~M 


FFT -'[U. 


C. Three wires 


HIGH PRESSURE 



FIG. 2. Schematic of rotating wheel apparatus. 


T <r=]vp FFT~'[r s „]- 


( 21 ) 


Consider three thermocouple wires of unequal diameters 
D\<D 2 <D 3 as shown in Fig. 1. The conservation expres- 
sion for each wire becomes 

dT \ 

dt ~ <0 i(T g ~T i), (15) 

dT 2 

-jf = <o 2 (T g -T 2 ), (1 6 ) 


dT 3 

-^- = a> 3 (T g -T 3 ), ( 17 ) 


with the additional constraint 


a — 


0)\ 


co l (o 3 \D,Z) 3 


D i 


m-2 


(18) 


Taking the FFT of Eqs. (15)— (17), squaring Eq. (16) and 
dividing by the product of Eqs. (15) and (17), one obtains 




gv f 2v ) 2 


TJ 3 j (f gv -f lv )(f gv -f 3v ) 


(19) 


Substituting r—T 2 /T 3 Ti an< ! solving for T gv , one obtains 


T = 

gv 


-b- yjb 2 -4ac 


2 a 


( 20 ) 


where 


a = r— a, 

b = 2aT 2v ~r(f lv + f 3v ), 
c = rf lv T 3v ~aT 2 2v , 

and r, e =FFT(|V|"7’ ? ). f u ^¥VY{\V\ m T { ). The inverse 
transform of T gv provides the reconstructed gas temperature 


III. EXPERIMENT 
A. Apparatus 

In the present experiment, thermocouple sensors are ex- 
posed to a reversing velocity air stream (<25 m/s) of varying 
temperature. In particular, the dynamic response of the ther- 
mocouple is measured for a periodic temperature profile of 
varying frequency. A rotating wheel configuration is used to 
deliver the test air stream to the proposed sensors. A similar 
experimental apparatus was described in detail by Elmore 
et al? and Forney et al. s 

A schematic of the rotating wheel apparatus used in the 
present experiment is shown in Fig. 2. As the wheel rotates, 
holes pass the two air supply tubes (3/4 in. i.d. copper) that 
allow slugs of hot (~55 °C) and cold (~30 °C) air to alter- 
nately enter the tube assembly mounted directly above the 
rotating wheel. The slugs of hot and cold air form a reversing 
air stream providing a periodic temperature profile covering 
a range of gas temperature frequencies from roughly 1 to 60 
Hz. 

The analog temperature signal is digitized with a Data 
Translation DT-2801 A/D board mounted in an expansion 
slot of an IBM AT compatible computer as shown in Fig. 2. 
ASYST software was loaded onto the hard disk of the per- 
sonal computer and this provided a flexible system for data 
storage, manipulation, and display. The true temperature pro- 
file of the airstream is measured with a constant current an- 
emometer (TSI 1054- A) and sensor (1210-T1.5). 

The three type K thermocouple wires were of diameter 
£>, = 50.8 fim (2 mil), D 2 = 101.6 jum, and D 3 = 203.2 
ixm, as shown in Fig. 1. The thermocouple junctions were 
fabricated by the Research Instrumentation Branch at the 
NASA Lewis Center. The wires were cut with a razor blade 
to produce a flat edge perpendicular to the wire axis. The 
wire segments were then mounted in a fixture so that the two 
faces of the junction were held together by springs. The junc- 
tion was then formed by laser heating. The laser used was a 
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FIG. 3. Velocity profile for reversing flow. Frequency of reversing flow 
/= 1 0 Hz. 



FIG. 5. Amplitude ratio of two thermocouple output-to-true gas tempera- 
ture. Nominal value for diameter ratio /3 0 =2.83 corresponding to manufac- 
turers’ diameters of D x = 2 and D 2 =4 mils. Optimum /3=2.55. 


KORAD model KWD Nd: YAG laser, operating at a wave- 
length of 1060 nm. The power settings depend on the wire 
size and material, but for the thermocouple described in this 
report, a pulse duration of 4 ms was used, delivering a total 
energy of approximately 2 J. A photograph of a junction 
formed in this manner is shown in an earlier article. 8 

B. Procedure 

ASYST software was developed to acquire temperature 
data sequentially from the thermocouple and constant current 
anemometer. The data are digitized for four channels at a 
sampling rate of 512 Hz per channel for a total sample time 
of 0.5 s. After data acquisition, the FFT is taken for each 
channel. The inverse transform of f gv by the ASYST code 
provided a reconstruction of the true gas temperature. 

IV. RESULTS AND DISCUSSION 

The velocity profile of the airstream was measured with 
a constant temperature anemometer. Figure 3 illustrates the 
velocity profile for a wheel frequency of 10 Hz at constant 



FIG. 4. Gas temperature profile compared to reconstructed temperature de- 
rived from a three-wire thermocouple. Frequency of temperature and flow 
reversal /= 2 1.5 Hz (w=135). Diameter ratio=1.0. The true gas tempera- 
ture profile is shifted to the right for clarity. 


temperature (heating tape off). For the positive velocities in 
Fig. 3, the gas reached a temperature of roughly 25-30 °C in 
the experiments below, corresponding to the cold slugs, 
while the negative velocity components were hot air slugs of 
approximately 50 °C. 

In the present study, the gas velocity and true tempera- 
ture were not measured simultaneously. Rather, it was as- 
sumed that the temperature measured with the anemometer 
was proportional to the gas velocity in the reversing flow 
field or V{t) a T g (t). In this case, \T g (t)\ m was substituted 
for | V(t)\ m throughout the analysis, where T g (t ) is the true 
gas temperature. The latter assumption was possible since 
V{t) appears in both the numerator and the denominator of 
Eq. (21). Thus, the amplitude of the gas velocity was not 
important while the phase angle of V(t) and T g (t) are iden- 
tical. 

A typical data profile constructed from Eq. (21) is illus- 
trated in Fig. 4. The uniform profile on the right is the tem- 
perature profile measured with the anemometer. The tem- 
perature profile on the left was reconstructed with Eq. (21) 
assuming that the exponent m = 0.5 in Eqs. (7) and (21) and 
a value of the diameter ratio a= 1 .0 for the three- wire ther- 
mocouple. The phase shift provided clarity between both sig- 
nals. 



FIG. 6. Amplitude ratio of two-thermocouple output vs gas temperature 
frequency of small wire based on maximum gas velocity. Angular frequency 
g)i =59 s _l . 
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FIG. 7. Amplitude ratio of three-thermocouple output-to-true gas tempera- 
ture. Nominal and optimum value for diameter ratio a 0 =1.0 corresponding 
to manufacturers’ diameters of D x = 2, D 2 = 4 and Z> 3 = 8 mils. 

The value of the exponent m — G.5 was chosen for the 
Reynolds number that appears in the Nusselt number corre- 
lation Eq. (5) and the wire natural frequency a ) n given by Eq. 
(7). At the maximum gas velocity of 25 m/s, the maximum 
Reynolds number for the three wires in the thermocouple 
assembly are 80, 160, and 320, corresponding to the smallest 
to largest diameter wire, respectively. Since the recom- 
mended values of the exponent m are 0.5 for Re>40 and 0.4 
for l=£Ress40, 6 we chose a constant value m = 0.5 in the 
present study. 

A. Two wires 

The amplitude ratio of thermocouple-to-gas temperature 
was determined using a two-wire combination with 
£>, = 50.8 (Am and D 2 = 101.6 /xm from the schematic in 
Fig. 1. The gas temperature was reconstructed with Eq. (13) 
assuming m = 0.5. The amplitude ratio was determined by 
averaging peak temperatures over six or eight cycles. The 
results are indicated in Fig. 5 as a function of the parameter 
P = (ti x l (02 = (£> 1 /Z> 2 ) 3/2 . In Fig. 5 the natural frequency 
<Wj=59 s is defined as the angular frequency at the maxi- 
mum gas velocity for the small wire D, = 50.8 ( an . It is 
clear that the amplitude ratio is strongly affected by (3 and 
that the optimum value is >8=2.55. 

It should be noted that the optimum (3=2.55 differs from 
an expected nominal value of £ 0 =2.83, where the latter is 
computed from the manufacturers’ wire diameters. These dif- 
ferences in (3 can be attributed to the empirical value of the 
exponent m = 0.5 chosen or to changes in wire diameter at 
the thermocouple junction because of the laser fabrication 
process. 

Figure 6 illustrates amplitude attenuation as a function of 
gas temperature frequency for fixed £=2.55. The two-wire 



combination requires no compensation for frequencies 
a»/wj<2.3 corresponding to frequencies /< 2 2 Hz in the re- 
versing flow field since 27 r/j = &> i = 59 s _I . 


B. Three wires 

The amplitude ratio of thermocouple-to-gas temperature 
was determined for the three-wire combination as shown in 
the schematic of Fig. 1. The gas temperature was recon- 
structed with Eqs. (20) and (21) with the exponent m = 0.5. 
The results are indicated in Fig. 7 as a function of the pa- 
rameter oi = ( o) 2 / ( 0 \ W3) = (D 2 / Z9 1 03) 3/2 . It is clear that the 
amplitude is strongly affected by a and that the optimum 
value is a^T.O. The latter value also corresponds to the 
nominal value calculated with the manufacturers’ listed wire 
diameters, as shown in Table I. 

Figure 8 illustrates amplitude attenuation as a function of 
gas temperature frequency for fixed a=1.0. The three-wire 
combination requires no compensation for frequencies 
u>ho\< 3.5 corresponding to frequencies /< 3 4 Hz in the re- 
versing flow field. 
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